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Summary 

A computing method for obtaining c.m.f. and average ground conductivity 
values from sets of ground-wave measurements taken on radials from an m.f. aerial 
should be valuable in assessing the performance of new arrays in the future. A suitable 
optimisation technique which can be employed for this purpose and could be realised as a 
computer program is described. 
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1. Introduction 

A m I e- of- thumb method has been in general use for 
estimating the ground-wave c.m.f.* and mean ground con- 
ductivity from sets of radial field measurements within a 
distance of about 5 km from an m.f. aerial. If £" is the 
field measured at distance d, the method is to plot log {Ed) 
against d on graph paper and to obtain the required results 
from the best-fit straight line drawn through the plotted 
points. This procedure can lead to very poor estimates of 
c.m.f., especially when the effective ground conductivity is 
low and the frequency is at the top of the band. The 
purpose of this report is to describe an improved method 
suitable for the assessment of aerial performance. It uses 
a non-iterative optimisation technique which can readily be 
realised as a computer program. This method can use any 
known theoretical attanuation function and optimises a 
scaling factor (equivalent to c.m.f.) and the ground con- 
ductivity value to give the minimum value of accumulated 
squared difference from the measured values of field 
strength. The final values of these quantities are obtained 
by interpolation. 
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it is assumed that fields will be measured over the 
ground (not at sea) and therefore it is possible to limit the 
expected values of conductivity and dielectric constant. 



(a) 



2. The attenuation function 

The ground-wave field strength for short distances (but 
greater than about 5X) can be expressed approximately as 



E=—^ Alp) 
d 



(VnT^) 



where A(p), the attenuation function. 



2 + 0-3p 
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In this approximation p is a numerical distance and i is a 
parameter. Both are defined below. This empirical 
function, due to Norton ' gives a good approximation to 
a more exact function as can be seen by reference to Fig. 
1 (a). A more exact expression derived from the theory of 
Sommerfield is 

A(5) =11-1- YJtFs . e~^ . erfc(-K/s)| 
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where erfc(z) is the complex error function complement Fig. 1 



which is tabulated,"^ and,S' = pe 



\b 



c.m.f. = cymo-motive force. For the purposes of this report it 
can be identified with the constant C in Equation (3). 



Theoretical curve fitting to measured results 

ia) Ed (Volts) as a function of distance 

- Approximate formula More exact formula 

ib) E {mV/m" ) as a function of distance 



Approximate formula 
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We see that the product £'d is equal to 7E^ atp = 0, 
(d = 0). 7E^ can therefore be identified with the aerial 
c.m.f. (Volts). Given the ground conductivity a (milli 
Siemen m~^ ), permittivity e (farad m~^) and the frequency 
/(Hz) the following quantities are required; 

x= (1-8 X 10')ff// (farad ?w~M 



tan b = 



(e+e„) 



p = — .cos 6. d. = kd 
xX 



The above formulae refer to vertical polarisation only. 
Horizontal polarisation is not covered in this report although 
a similar set of procedures can be devised for it. 



3. Ground permittivity 

The permittivity of different types of ground is, of 
course, not a well defined function of conductivity. In 
temperate climates, however, there is a tendency for per- 
mittivity to increase with conductivity which is probably 
due to increasing water content in the soil. In cities and 
built-up areas 'real' conductivities and permittivities must 
be replaced by 'effective' values both of which are low. 
Rather than use a fixed typical or mean value of per- 
mittivity (as is done in CCIR Recommendation 368-1 for 
example), improved accuracy can be secured by expressing 
the relative permittivity as a single valued function of con- 
ductivity, Thus 



20(1 -e' 



-0-23CT1 



(2) 



a = conductivity in mS, m 

This formula will not necessarily give the correct value of 
permittivity in a place of a fixed value but its justification 
is that on the average smaller errors will be incurred if it is 
adopted. Fig. 2 shows a plot of this function together 
with some typical examples of points corresponding to 
various types of terrain. 



4. IVIeasured£'c? values 

Restating the basic formula in a slightly different form: 

V = Ed=CA{kd) (Volts) (3) 

It will be assumed that a set of measured values oi Ed is 
available which has been taken on a radial from a vertical 
transmitting aerial within a range of approximately 5 to 
200 wavelengths. The given data are therefore a set of 
points {d,^ Vj.) ;■ = 1, 2 . . . n at the frequency of the test 
transmission. Owing to practical measuring difficulties the 
number of measurement points may be limited to about a 
dozen and often less. 
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Fig. 2 - Assumed relationship between conductivity and 
relative permittivity for overland paths 

5. Optimisation method 

5.1. Step 1 

Chose \os\, a set of values of conductivity in the 
range 0-1<o<20, say. It is convenient, as far as possible, 
to have roughly equal increments between the chosen values 
of a and a set of about fifteen to twenty is sufficient. 
Calculate corresponding values of k giving the set \kX 
From thesevalues a set of values of the attenuation function 
can be found from the formulae in Section 2. 



\/K,,\^[Mk^d,)\ 



These values and the points (d^, V^) form the arrays on 
which the optimisation procedure is based. 

5.2. Step 2 

\cA is next derived where 



C. 



E 



r=1 



Asr^^r/ 



'E'^^^'^ 



r=1 



This process scales the theoretical F-curves for each of the 
conductivity values by adjusting the multiplying constant 
C in expression (3) to minimise the sum of squared dif- 
ferences. 

5.3. Step 3 

This step involves the calculation of the sum of 
squared differences between measured and calculated values 
of F. 
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Thus 



3iving|5^|. 






5^= > ^[F,-CArI 
r=1 



7.2. Step 5 

Given the optimal value of a, obtain the optimal 
value of C by linear interpolation or by the application of 
Step 2 for the single optimal value of a. 

8. Example 



6. Functional dependence on conductivity 

The multiplier C in Equation (3) is to be thought of as 
equivalent to the c.m.f. of the aerial system which is being 
measured. In reality this multiplier has a fixed value 
determined by the aerial system but in the optimisation 
process its value changes as different values of conductivity 
are assumed. In the present context it may therefore be 
considered as dependent upon the conductivity. It follows 
that both the multiplier C and the squared difference 
variable are functions of the conductivity only which means 
that a one-dimensional interpolation may be used to 
obtain the best estimate of conductivity and c.m.f. as 
explained in the next section. 



Fig. 1 shows the result of applying the process to six 
radial field strength measurements taken from a service 
broadcast aerial operating at a frequency of 1457 kHz. 
From the measurements the aerial c.m.f. (in the radial 
direction) is found to be 221-9 volts and the effective 



ground conductivity to be 3-07 mi' m' 



This is a fairly 



low value of conductivity but the path of the radial is 
known to pass over undulating ground which is well 
covered with banks of trees and buildings and this no 
doubt accounts for this fact at the relatively high frequency. 

The dotted curve in Fig. 1(a) shows a plot of the 
more exact function mentioned in Section 2 derived from 
curves given in Reference 1. 



7. Interpolation 

The value of conductivity which corresponds to the 
minimum of 6^ can be found to good accuracy by means of 
parabolic interpolation. The interpolation involves the 
point (5|^ CT|) giving the minimal calculated value and the 
neighbouring points on either side. It can be shown that 
the interpolated value of a will be an accurate estimate of 
that value for which 5^ is a minimum. Having found this 
optimal value of a, the best estimate forC (the c.m.f.) may 
be obtained by linear interpolation or by applying Step 2 
(Section 5.2) for this single value of conductivity. 



9. C.M.F. ground-wave patterns of aerials 

One of the principal applications for the process des- 
cribed in the preceding sections is the calculation of c.m.f. 
ground-wave patterns of vertical m.f. aerials from a series of 
measurements. This application will be particularly im- 
portant when new directional arrays are being set to work 
and the amplitude and phases of their drive currents are 
being critically adjusted. The advantage of the process is 
that it minimises ground attenuation effects in the estima- 
tion of aerial performance. 



7.1. Step 4 



M 



in |5j> is found and, if this is a single member of 
the set, 5| say, the points (5f_i, 0\_-^) {8f , ffj) and 
(5|^+l, CT|+i) are noted and interpolation carried out. The 
optimal value of a is given by: 



2a = 



iaf-al,){8l,-8f)-{al, 



■a^){S{ 



2 
1+1 ' 



•5^ 



io. 



1 + 1 



) (5|_i -5|')-(a| 



1-1 



a,) (5 



•5f+i) 



If there are two sequential members of \83 h 

:hosen aso?. 



having 
the same minimal value either may be chosen as5| . The 
likelihood of more than two sequential members having the 
same value is negligible. If the minimal value belongs to 
either end member (but not both) it can be taken as the 
actual minimum. In this event, however, note should be 
made of the fact that the range of conductivity values may 
be inadequate. If nominally equal minimal values are 
possessed by separated members it should be inferred that 
the measured results are too erratic for smoothing to be 
applied. 

* The possibility of this occurrence will clearly depend on the 
number of significant figures which are available. 



10. Ground conductivity 

In the author's opinion calculation of m.f. ground-wave 
services often assume values of effective ground conduc- 
tivity which are too high when the operating frequency is at 
the top of the m.f. band (above 1 MHz). Ground undula- 
tions, trees, building and other irregularities cause increased 
attenuation whereas predictions are often based on surface 
stratum values of conductivity and measurements made at 
low frequencies. Application of the optimisation process 
to measured fields from aerials could lead to improved 
estimates for the average conductivity to be assumed in 
future work. 



n. Conclusions 

An explicit non-iterative optimisation process has been 
described for obtaining the ground-wave c.m.f. values and 
local effective ground conductivity of vertical m.f. aerials 
from measurements. The accuracy of the method as a 
curve fitting process is believed to be good and its applica- 
tion in assessing the performance of new m.f. arrays in the 
future will, it is hoped, be valuable. 
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